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1 Introduction and Formulation of Results 



O 

■ Abstract 

We explicitly calculate the Janossy densities for a special class of finite determinantal random 
point processes with several types of particles introduced by Prahofer and Spohn and, in the 
full generality, by Johansson in connection with the analysis of polynuclear growth models. The 
results of this paper generalize the theorem we proved earlier with Borodin about the Janossy 
' densities in biorthogonal ensembles. In particular our results can be applied to ensembles of 

random matrices coupled in a chain which provide a very important example of determinantal 
ensembles we study. 

> 

On 

o 

1.1. Janossy densities. The main goal of this paper is study the Janossy densities in a special 
class of finite determinantal random point processes introduced recently by by Prahofer and Spohn 
(|32j). and, in the full generality, by Johansson ( in connection with the analysis of a certain 
class of polynuclear growth models (the papers and ^H] considered sets of non-intersecting line 
Q^' ensembles with fixed initial and final points). A similar framework was introduced by Okounkov 

and Reshetikhin an d Ferrari and Spohn to analyze the 3D Young tableaux via non- 

intersecting line ensembles. The distribution of the eigenvalues of random matrices coupled in a 
chain (see |0j, [22], |Sj P] also falls into this class. The preprint (T2J shows that the 

the statistics of the random field built from the PNG model is related to the edge properties of a 
Gaussian multi-matrix model. 

The term Janossy densities in the theory of random point processes was introduced by Srinvasan 
in 1969 ([HE]) who referred to the 1950 paper by Janossy (|18j) on particle showers. The classical 
reference nowdays is (jTj, see chapters 5 and 7). We postpone till section 1.2 the formal definition 
of the Janossy densities for the determinantal ensemble we are interested in. In order to give an 
uninitiated reader a glimpse of what is going on let us consider a random point process with one 
class of particles on R (i.e. a probability measure on a space of locally finite point configuration 
in R) and assume that all point correlation functions exist and locally integrable (for all practical 
purposes the reader can think about a Poisson random point process for now). For a finite interval 
/ C R we can define the k-point Janossy density Jk,i{x\, . . . , Xk) for x±, . . . , xj. E I as 
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<Jk,l( x ii ■ • ■ i x k) dxi = Pr{ there are exactly k particles in / and there is a particle in each 

of the k infinitesimal intervals (ccj, x; t + dxi), i = 1, . . . , k}. (1) 
A very useful property of the Janossy densities is that 

Prithere are exactly k particles in 1} = — [ J k i(x\, . . . , x k )dx\ ■ ■ ■ dx k . (2) 

In particular if we know that the number of the particles in a configuration is finite (which is 
often the case in Random Matrix Theory and other applications) one can use Janossy densities to 
study the distribution of the largest /smallest particles (eigenvalues). Let us order the particles in a 
finite configuration Ai > A2 > A3 . . . > A n . Then the probability density of the k— th largest particle 
can be expressed in an integral form as 

Pr(A fc G (s, s + ds)) = ( 1 / Jk,(s,+oo)(.xi, • • • > %k-i, s)dx\ ■ ■ ■ dx k _i ] ds, (3) 

\(k-l)\ J {Si+oo) )k-i J 

where in @ we put x^ = s. The equivalent form of © is 
Pr(A fe > s) - Pr(A fc+ i > s) = Pr(#[s, +00) = k) = — [ Jk,{s,+oo){xi, • • • , x k )dxi ■ ■ ■ dx k . (4) 

J(s,+oo) k 

Janossy densities are equally useful if one is interested to study the local statistical properties of 
the distribution of particles in the bulk of the configuration. 

In ([I]) we studied the Janossy densities in a so-called biorthogonal ensemble (see e.g. [S] and 
|40| ) which can be described as a random configuration of n particles in a one-particle space X (for 
simplicity we will consider the case X = K) with the joint probability density (with respect to the 
Lebesgue measure in the case X = M) 

p n (xi,...,x n ) = — det(/ i (a; J O)?j=idet(0i(a;j))?j=i, (5) 

where {fi(x), (j>i(x), i = 1, . . . ,n} are some complex-valued functions on M. Such ensembles were 
extensively studied in random matrix theory ([H], HZj), directed percolation and tiling models ([2P], 
|21j . P2J, models of uniform spanning trees and forests on graphs (jSJ, |25| ) and representation 
theory (0), among others. For the biorthogonal ensemble one can calculate explicitely the k— point 
correlation functions 

p[ n \x 1 ,x 2 , ...,x k )= n ' [ p n (x 1 ,...,x n )dx k+1 dx k+2 ---dx n . (6) 
(n - k)\ J R n-k 

It is a standrad fact that k— point correlation functions have a determinantal form in this case 

p k n) (xi,x 2 , ...,x k ) = det (K(xi,Xj)) ij=htk , (7) 

n 

K(x,y) = J2hx)My), (8) 
i=l 
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where {/j, (pi, i = 1, . . . , re} are biorthogonal bases in Span{fi, 1 < i < n} and Span{(pi, 1 < i < n}, 
i.e. 

J fi(x)<j)j(x)dx = 5ij. (9) 

Then the Janossy densities also have the determinantal form (see [J], p. 140 or [2], Section 2) with 
a kernel C 1 : 

J k j(x 1 , . . . ,x k ) = coresi(/)det(£ 7 (xi, 2^)^=1,...^), (10) 

where 

L 1 = Ki{ld - KiY x , const(I) = det(Id-K I ), (11) 

and the kernel of the integral operator Kj is the restriction of the kernel of K to / x /, i.e. Kj(x, y) = 
Xi(x)K(x, y)xi{y)i where xi 1S the indicator function of /. The probabilistic interpretation of cost(I) 
is that of the probability to have no particles in /. 

The formula for the Janossy kernel C 1 for such ensembles was explicitely calculated in (@[). It 
was proved that the kernel C 1 can be constructed according to the following rule: 

1) consider the ensemble with one-particle space X replaced by X \ I (i.e. the density of the 
distribution is still given by the same formula, only now it is defined on {X\I) n , rather than on X n ; 
naturally the normalization constant changes). 

2) calculate the correlation kernel on {X \ I) x [X \ I) using Q), (i.e. the pairing to be 
considered on the functions on X\I). 

3) extend the correlation kernel to I x I (since for M = 1 the correlation kernel is expressed in 
terms of fi, (ftj there is no ambiguity in how to extend it to I x I). 

In the special case of a polynomial ensemble of Hermitian random matrices (which corresponds 
to fj(x) = <pj{x) = xi~ l exp (— ^y(x)) , one can see that the Janossy kernel can be written as the 
Christoffel-Darboux kernel C J (x,y) = X]j=o Pj( x )Pj(y) built from the orthonormal polynomials with 
respect to the weight fi(dx) = exp(—V(x))xi c (x)dx. As one can see © and (HJ) lead to a nice limiting 
expression for the distribution of the (appropriately rescaled) k— th largest eigenvalue in the limit 
re — > oo provided one can handle the asymptotics of the orthogonal polynomials with respect to the 
weight of the form e~ nV ^ X(-oo,s)( x ) (where s = s n is at the "right edge of the spectrum"). Recent 
remarkable advances in the application of Riemann-Hilbert problem technique to the asymptotics of 
orthogonal polynomials (see e.g. [2]) suggest to us that this might be possible to. In particular in the 
GUE case the problem is thus reduced to the calculation of the asymptotics of the re— th orthonormal 
polynomial with respect to the weight exp(— £ 2 )X(-oo,s n )( x )i where s n = 2 1 / 2 n 1/ ' 2 + x2 -1 / 2 n~ 1 / 6 and 
x is a fixed real number. The limiting formulas suggested by © and @ would be, in our opinion, 
simpler than the ones currently used that require k differentiations of the Fredholm determinant 
of the infinite-dimensional integral operator det(l + zK/ St+oa \) with respect to the parameter z. In 
a joint paper with Borodin ([I]) we were able to calculate explicitely the asymptotics of the k— th 
smallest eigenvalue in the standard Laguerre (Wishart) ensembles (see also [42 ). 

In (jSHJ) we proved that the same recipe (with simple alterations) applies to pfaffian ensembles 
given by the formula ([33]) [3T)j . see also [IT] ) 

p(yi,---,V2n) = v-det(/ij(y fc ))i.fc=i,.-,2n. P/(e(%>yfc))j,fc=i,...,2n, (12) 

where hi, ... , hm are complex-valued functions on a measure space (Y, d\(y)), and e(y, z) is a skew- 
symmetric kernel, e(y,z) = —e(z,y). For the definition of the pfaffian of a 2n x 2re skew-symmetric 
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matrix we refer the reader to ( [IB] ) . The pfaffian ensemble (|12|) contains the biorthogonal ensemble as 
a specail case. Some other interesting examples of the pfaffian ensembles include (3 = 1 (orthogonal) 
and (5 = 4 (symplectic) polynomial ensembles of random matrices (see |36j ) as well as the examples 
appearing in works on random growth models and vicious random walks (we refer to ^H] and |29j 
and the references in there). The name pfaffian for the ensemble 1)12(1 comes from the fact that both 
point correlation functions and Janossy densities have the pfaffian form (see [22], j30j ) . 

1.2. Determinantal ensembles with several classes of particles. Now we are ready to turn to 
the class of determinantal ensembles introduced in the papers by Johansson (^HJ) an d Prahofer and 
Spohn (|32j)- Let (X, fj,) be a measure space, fx, f'2, ■ ■ ■ f n , 4>i> 4>2, ■ ■ ■ , 4>n - complex- valued bounded 
integrable functions on X, and 31,2 y),92,3(x, y), . . . , gM-i,At(%> y) - complex- valued bounded inte- 
grable functions on X 2 = X x X with respect to the product measure //® 2 = fi x \i (in principle the 
above assumptions on fj, gii+i, (pi, i,j = l,...,n, I = 1, . . . ,M can be weakened). Suppose that 

D nr(r {1) tW-T (2) ■ ' T (M) T^l 

t>n,M V-^i ? • • • > •»<n ) -^l 1 ■ ■ ■ 1 ""n ' " ' ' ' 1 ' ' ' ' ' n I 

M-l 

— det(/i(x«))^. =1 J] det(g ljl+1 (x?\ x« + x )))^ =1 det(^-(xf )))^ (13) 



M-l 

, det(/i(x«))^. =1 



defines the density of a M x n- dimensional probability distribution on X = X x • • • x X with 
respect to the product measure n® Mn . One can view the configuration 

x = (x^, . . . , x^; x^\ . . . , xffl; x[ M \ . . . , ii M ') as the union of M configurations, namely the 
first floor configuration = (x^\ . . . ,x&), the second floor configuration x( 2 ) = {xf \ . . . ,Xn^), 
etc. In other words we can call the particles of the first floor configuration - particles of the first class, 
the particles of the second floor configuration - particles of the second class, etc. In the original papers 
(.32 ), (^Hl) the formulas of the type (fT3*|) appeared as Prahofer, Spohn and Johansson considered 
sets of non-intersecting line ensembles with fixed initial and final points. 
The normalization constant in 1)13(1 (usually called the partition function) 

„ M-l M n 

Zn,M= / m det(/ 4 (xj 1) ))^ =1 J] det( 5/ , /+1 (xf\x( /+1) ))^ =1 det(^(xf ) ))^ =i nn^( d2; f ) ) 
JxMn i=i i=i i=i 

(14) 

can be shown to be equal (n\) M det(^4), where the n x n matrix A = (A//c).,',fc=i,...,n is defined as 

. M-l M 

A jk = / fj(x {1) ) J] 9i,i + i(x (l) ,x^)Mx (M) ) II Kdx^). (15) 

J XM 1=1 m=l 

We assume that the matrix A is invertible. 

It is easy to see that the biorthogonal ensemble (J5J) is a special case of (|T3*)l corresponding to 
M = 1. What is more interesting is that the ensemble ()13l) in the case of two classes of particles (i.e. 
M = 2) is a special case of the pfaffain ensemble 1(12(1 discussed at the end of section 1.1. Indeed let Y 
be the disjoint union of two identical copies of X, Y = X\ \_\ Xi, and the restriction of the measure 
A on each copy of X given by fi. Suppose that for 1 < i < n the restriction of hi on X\ is given 
by fi and the restriction of hi on X2 is identically zero. Similarly, suppose that the restriction of 
h n +i on X\ is identically zero, and the restriction of h n+ i on X2 is given by <fii, i = 1, . . . ,n. Finally, 
suppose that the kernel e is identically zero on X\ x X\ and X2 x X%, and e(xi, X2) = —e(x2, x\) for 
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x\ G X\, %2 G -X2. Let us define a kernel g on X x X which takes the same values as e on X\ x X 2 . 
Then the formula (|12|) specializes into (|13j) . M = 2 . 

For the ensemble (fT3"|) one can explicitly calculate (k%, k 2 , ■ ■ ■ , fc/wO-point correlation functions 

„ M n 

Pfc 1 ,...,fc M (xi 1) ,...,x[, 1) ;...;xi A//) ,...,x^ ) ) := / Pu,m{x) TT(n!/(n - fc z )!) TT dn(xf), (16) 

JX 1=1 j=k+i 

where k = &i H 1- /cjvf , < kj < n, and show that they have the determinantal form (|19j. |32| . see 

also jH]) 



«i,...,*;A/ri j • • • j }z\ ' * * * ' 1 •'■■■•>- h k M ) 
= Aet{K n ' M {l,xf\m,x^)) hm= x^ M , l^k^i^m^*™- (17) 
To define the kernel /C we introduce the following notations for the convolutions: 



9i,l+i*9i+i,i+2{x,y) ■= gi ) i +1 (x,z)gi + x > i +2 (z,y)dn(z) (18) 

Jx 

9l,m ■= 9i,i+i * • • • * 9m-l,m, I < I < m < M (19) 

9l,m := 0, 1 < m < / < M. (20) 

We will use similar notations for the integrals J x fj(x)gi tm (x,y)d^(y) and J x g m ,M (x, y)4>s(y)d/J<(y) , 
namely 



(fj*9l,m)(y)= fj(x)gi, m (x,y)dfi(x), j = l,...,n, (21) 
Jx 

(g m ,M * 4>s)(x) = / 9m,Ai(x,y)(ps(y)dn(y), s = l,...,n. (22) 



The kernel /C n,M : ({1,2,... , M}) x X) 2 1— > C is a M x M matrix kernel given by the following 
expression 

n 

lC n > M (l,x;m,y) = -g lm (x,y) + ^ {g hM * ^{x^A'^tfj * g 1>m ) (23) 

i,i=l 

(to simplify the formulas we adopted the convention f, L * <7i i7n = /j for m = 1 and * </>j = 
for I = M). Usually we omit the dependence on n and M in the notation of the kernel if it does not 
lead to ambiguity. 
Remark 

Repeated applications of the Heine identity 

— j det(ipi(xj))ij = i t ..., n det(ipi(xj))i t j=i,..., n TT dfx{x k ) = det( / <J>i{y)4>j(y)dn(y)Y ,, 
n! JX" ;~=i Jx ,J 



to implies that the joint distribution of the h, ■ ■ ■ , lm- floors configurations x( ll \ . . . , x( lm \ 1 < 
h < ...l m < M, is again 0/ determinantal form f73J) with fa = fi * gi^, g\ )2 = gi u i 2 , £12,3 = 
9h,h,--- ,9m-i,m = 9l m -i,l m i h = 9l m ,M * <f>j, M = m and A = A. 



5 



If X C M and fi is absolutely continuous with respect to the Lebesgue measure, then the proba- 
bilistic meaning of the (k\, . . . , &m) -point correlation functions is that of the density of probability 
to find a particle of the first class in each infinitesimal interval around points x^\ . . . , x^\ a particle 

(2) (2) 

of the second class in each infinitesimal interval around points x\ , . . . , x^ , etc. In other words 

pki,...,k M ( x ? i ■ ■ ■ >4?; • • • ;x i M) ' • • ''CVrf' ' ' ' ^ dx £) = 

Pr{ for each 1 < I < M there is a particle of the Z-th class in each of the ki intervals 
{xf,x ( P+dxf), l<i<h). 

On the other hand, if \x is supported by a discrete set of points, then 
Pfci,...,fc M ( x i ) • • • ; x ki'i ■ ■ ■ > x l i • • • i x kM^( Xl ) * " ^ X kM^ = 

Pr { for each 1 < I < M there is a particle of the Z-th class at each of the sites xf\ i = 1, . . . , ki}. 

In general, random point processes with the point correlation functions of the determinantal form 
are called determinantal (a.k.a. fermion) random point processes (|37|). 

The Janossy densities Jit 1 ,/ i; it 2 ,/ 2 ;...,ifc M ,/ M (a;^, . . . , x£; . . . ; x[ M \ x^), < k < n, I = 
1, . . . , M, k\ + . . . + ku = k < M x n, describe the joint distribution of the first class particles 
in I\ , second class particles in I2, . . . , M-th class particles in 1^ , where I\ , I2, ■ ■ ■ , Ik are measurable 
subsets of X. For the ensembles with a finite number of particles Janossy densities can be obtained 
from the joint probability distribution of the particles by integration, in particular for (|13|) we have 



Jki,I 1 ;k 2 ,l2\---,k M ,I M {x 1 , . . . , X k \ . . . ) X± ' 



(1) Jl) .JM) JM), 

M n 

PnM x )I\( nl /^-kV.) n Mx?)> ( 24 ) 

(X\h)*ix-x(X\I M ) k M i= ^ A +1 

One can say that the Janossy density Jk 1 ,h;k2,l2;-;k M ,lM( x i ; • • • > x ki'' ■ • • ! x i M) > • • • > x k^) § ives the 
joint density of the distribution of k\ first class particles in I\, k<i second class particles in I2, 
ku M-th class particles in Im (under the assumption that there are no other particles of the first 
class in I\, no other particles of the second kind in I2, etc). The Janossy densities differ from the 
the conditional probability densities by the normalization: the Janossy densities are normalized in 
such a way that the whole mass is not one, but rather 



If M kM 

V ... fcAf i L TkM <7kuh;-;k M ,iM( x i\- ■ ■ ,4?; • ■ ■ . x i M) > ■ ■ ■ > 4m } ) II II d ^ x n) 

1- M • J l 1 X...Xl M l=Hi=l 



M k M 

I T, r. ,. r (nP ^ 

hl---k M \ 

Pr { for each 1 < I < M there are exactly k\ particles of the I— th class in //} 



Let 4 , . . . , be some distinct points of I\ , x^ , . . . , scj^ - some distinct points of I2, ■ ■ ■ , x[ M ^ , . . . , x\^ 
some distinct points of Im- If X C R and \i is absolutely continuous with respect to the Lebesgue 
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measure, then 



Jkrjrtej2;...,k M j M ( x i i ■ ■ ■ >4i ; • • • ; 4 > • • • '4«V(^i } ) • • • ) = 

Pr { there are exactly fci particles of the first class in I±, . . . , ku particles of the M-th class in Im, 
so that for each 1 < I < M there is a particle of the Z-th class in each of the ki infinitesimal 
intervals (xf\ xp + dx^), i = 1, . . . , fc;}. 

Similarly, if /i is discrete, then 

7t r t r t r fx (1) -x (M) x (M) )ufx (1) V--ufx (M) ) - 

Jk 1 ,I 1 ;k 2 ,h;~;k M ,lM\ x l i ' ' ' > x fci ' ' ' ' ' x l ' ' ' ' » x fc M / ^ k M ) ~ 

Pr { there are exactly fci particles of the first class in I\, . . . , Jzm particles of the M-th class in Im, 
so that for each 1 < I < M there is a particle of the l-th class at each of the k\ sites 
xf\ i = l,...,k}. 

See [7], [I] and [HZ] for additional discussion. It is instructive to compare the probabilistic interpre- 
tation of the Janossy densities with the probabilistic interpretation of the correlation functions given 
above. For determinantal processes the Janossy densities also have the determinantal form (see jjj, 
p. 140 or [2], Section 2) with a kernel C 1 : 

% r u r (x (1) x (1) - x (M) x (M) ) 
= const(J)det(£ x (Z,x^ ) ;m,x^ ) ));, m =i,...,A/, i<i ; <fe ; , \<i m <k m , (25) 

where 

C 1 = Kj{Id-Kj)-\ (26) 

and the notations X, const{Z) and K,% are explained in the next paragraph. One can point out 
that we already have seen formulas of this type when we discussed a special case of biorthogonal 
ensembles. 

The integral operator 1C acts on a Hilbert space 7i, which is the orthogonal direct sum of M 
copies of L 2 (X, /i), i.e. TC = L 2 (X, /u) ... L 2 (X, fi). Let X be the disjoint union of M identical 
copies of X, in other words X = X± \_\ . . . \_\ Xm, where each Xi, 1 = 1,..., M, is a copy of X. One 
can think of X\, 1 < I < M, being the l-th floor in our particle space. Extending the measure 
fi in a natural way to X and denoting the extension by \im we can view 7i as the Hilbert space 
L 2 (X, //m)- For I\ C X, . . . ,Im C X, we construct a subset of the particle space X , denoted by 
X, in such a way that the intersection of I with X\ is equal to I\, 1 = 1,... ,M. Let us denote 
by K,j the restriction of the integral operator K to L 2 (T,hm) = L 2 {Ii) . . . L 2 (Im) (in other 
words we restrict the kernel /C to X x Z.) The normalization constant const{I\, . . . , Im) = const{Z) 
is given by the Fredholm determinant const(Z) = det(Id — ICj) of the operator ICj (for the definition 
of the Fredholm determinant we refer the reader to |H4] . j.S5j ) . The probabilistic meaning of the 
normalization constant const (Z) is that of the probability to have no first class particles in I\, no 
second class particles in Ii, no M-th class particles in Im- 

Remark 

Strictly speaking, to view the ensemble as a determinantal random point process we have to 
consider it as a distribution of M x n identical particles on a one-particle space X , so that with 
probability 1 there are exactly n particles in each X\, 1 = 1,..., M, and the joint distribution of the 
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particles in X\, X2, ■ ■ ■ , Xm is given by hl'J\) . Then the (ki,k2, ■ ■ ■ , &m) -point correlation function is 
nothings else but a usual k-point correlation function, k = k\ + ■ ■ ■ + kMi an d ki is just the number 
of arguments of the k-point correlation function that belong to Xi, I = 1, . . . , M. In particular one 
can easily see that 



£ / J C(t.« i ^ WK m,,)*ft f )-K(*^m,,) + (m-*VC(*, a r,m,.)+2(m-t) ai ,.(«,.) (27) 



M 

and an easy generalization of the Dyson-Mehta lemma (\2T$ , Theorem 5.2.1) claims that if we in- 
tegrate out ( over X ) a variable in the determinant of the k x k matrix with the correlation kernel 
K. we obtain (up to a trivial combinatorial coefficient) the determinant of the (k — 1) x {k — 1) ma- 
trix with the same kernel. To the reader familiar with the original Dyson-Mehta argument we note 
that the terms with the factor (m — k) vanish since for any permutation a € Sf. one (trivially) has 
Sj=l(°"(0 ~~ + 1)) = 0) where we put a{k + 1) := cr(l). 

1.3. Formulation of the main result. If one can prove that the analogue the recipe for the Janossy 
kernel described above in the biorthogonal case (after formulas applies to the ensemble 

(|13l) for general M, the Janossy kernel would have the following form: 

n 

£ J (Z,a;;m,y) = -gf m (x,y) + J2 *c cp i )(x)(A c )- j 1 (f j * c <? 1>m )(y), (28) 

£,j=l 

where 

9l,l+i *c9i+i,i+2(x,y) := / gi,i+i(x, z)gi +1 j +2 (z, y)dfi(z) (29) 

/ T c 

J h+i 

9l,m '■= 91,1+1 *c ••• *c 9m-l,m, I < I < m < M (30) 
g? m := 0, 1 < m < I < M, (31) 



(fj *c 9i, m )(y) = fj{x)gl m (x,y)dfi(x), j = l,...,n, (32) 
J If 

(9m,M *c(f>s)(x) = / g c mM {x,y)4> s {y)d^i{y), s = 1, ...,n, (33) 

/ T c 

M 

. M-1 M 

A% = I f 3 {x {1) ) II 9l,l + i(x«\x( l+ V)<t> k (x {M) ) II »(dx (m) ), (34) 

Jlfx...xI c M l=1 1 



m=l 



(the notation I c stands for the complement of a set /, we also remind the reader that we use the 
convention fj * c g\ m = fj for m = 1 and gf M * c (p s = <j) s for / = M.) 

Throughout the paper we assume that the matrices A (defined in (|15j) ). A c (defined in (|34|l ) and 
A 1 (defined below in ©) are invertible. The main result of this paper is 

Theorem 1.1 The Janossy kernel in the determinantal ensemble kl&\) is given by the formulas 
for all M. 
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To finish this section we introduce the notations for the convolutions over I\, . . . , Im- 

9i,i+i *i 9i+i,i+2{x,y) := / gi t i +1 (x,z)gi +l! i +2 (z,y)d(j,{z) (35) 

Jh+i 

(fj *i9i,m)(y) = / fj(x)gi, m (x,y)dfi(y), j = 1, . . . ,n, (36) 
Jl x 



(g m ,M *x<f>s){x) = / g m ,M( x iy) ( l ) s(y)dn(y), s = l,...,n, (37) 
Jim 

. M-l M 

a% = / n 5M+1 (*« ^ +i >)^(*( m )) n ^x^). (38) 

Jhx...xI M l=1 m=1 

The case M = 1 of the Theorem was proven in 0J. The case M = 2 ( "two-matrix model") 
follows from our proof for the pfaffian ensembles ()12|) given in ( ^36 ) . The case M = 3 ( "three-matrix 
model") will be proven in section 3 as a "warm up". The proof of the Theorm for general M is given 
in section 4. We devote section 2 to some additional examples of the determinantal ensembles (|13[). 

2 Examples of Determinantal Ensembles 

One Matrix Models. Unitary Ensembles 

Let M = 1. In the special cases X = R, /j(x) = 0j(x) = a^ -1 , and X = {C | |z| = l},fj(z) = 
4>j{z) = z J , such ensembles are well known in Random Matrix Theory as unitary ensembles, see 
|27| for details. An ensemble of the form © which is different from random matrix ensembles was 
studied in (2HJ - 

Matrices Coupled in a Chain 

Consider the chain of M complex Hermitian n x n matrices with the joint probability density 
(with respect to the M x n 2 -dimensional Lebesgue measure dA{) given by the formula (see [2*7] . 

M) 

F(Ai, ...,A M )= const(n,M)exp (-Tr^V^Ax) + V 2 (A 2 ) + ■■■ + V M -i(A M -i) + ^V M {A M )]\ 
x exp (Tr{ci^iA 2 + c 2 A 2 A 3 + ■■■ + c M -iA M -iA M }) . (39) 

The case M = 1 corresponds to the one matrix model discussed above. Let us denote by \f , , . . . , A„ 
the eigenvalues (all real) of A[, I = 1, . . . , n. The probability density of the joint distribution of the 
eigenvalues of A\,... ,Am with respect to the Lebesgue measure on M. Mn is given by the formula 

Pi*?,- ■ ■ , 4 M) ) = ^- exp {- £{ W } ) + V 2 {^) + . . . + VmM^) + \vm(4 M) V 



M-l 

x II - 4^ - ^ II det(exp( Q ,f)xf 1 )))^ =1 . (40) 

l<i<j<n 1=1 

Writing the Vandermonde products in (|40|) as determinants we arrive at the expression of the form 
(|13|) , In should be noted that while in the case of one- matrix polynomial ensembles the existing 
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Riemann-Hilbert problem machinery (see e.g. [H]) could, in principle, provide the desired asymptotics 
of the orthogonal polynomials with respect to the weights exp(— nV(x))xi% (x), this is not yet the 
case for multi-matrix models. 

Non-Intersecting Paths of a Markov Process 

We follow |24| . |23| . Let Pt,s(x,y) be the transition probability of a Markov process on K. 
with continuous trajectories and (£i(£), £2(1), ■ ■ ■ 5 Cn(*)) " n independent copies of the process. A 
beautiful classical result of Karlin and McGregor states that if n particles start at the positions 
a?j < x% < . . . < Xn \ then the probability density of their joint distribution at time t\ > 0, given 
that their paths have not intersected for all < t < t\, is equal to 

n tl (x?\ . . . = det^^.xf))^! 

provided the process £2(^)5 • • • >£n(i)) hi ^ n has a strong Markovian property. To understand 

the above formula better one can consider first the case of two particles and use a standard reflection 
trick to check that the result is correct. The most general combinatorial form of Karlin-McGregor 
theorem is known as the Gessel-Viennot theorem (HH), we refer the reader for the additional dis- 
cussion to jnn]), section 2.7. 

Let < t\ < t 2 < ■ ■ ■ < tM+i- The conditional probability density that the particles are in the 
positions x± < x% < . . . < xffi at time t\, at the positions x^ < x^ < . . . < i{, 2 ' at time t2,—, 
at the positions x^^ < x 2 M ^ < . . . < x^^ at time tM, given that at time £m+i they are at the 
positions x^ I+1 ^ < x^ I+1 ^ < . . . < x^ I+l ^ and their paths have not intersected, is then equal to 

1 M 



1=0 



where to = 0. One can easily see that (|41jl belongs to the class of ensembles (j!3j) . As an interesting 
related example we refer to the random walks on a discrete circle (see |13| and |19j . section 2.3) 

Finally we refer to ( |19j and |32j . specifically to the formulas (1.17)-(1.19), (3.15)-(3.16) in the first 
reference) for an example of a determinantal ensemble ()13|) appearing in the analysis of a polynuclear 
growth model. 

3 Case of Three Classes of Particles. 

We devote this section to the proof of our main result in the special case of three classes of particles. 
In the case M = 3 the formulas (|13H15j) . (|23j) have the following form: 



D J r {1) r W-r (2) T ( 2 )- r ( 3 ) T (3h 
fn,>3V'*'l 5***5 M n 5 -^l 5***5 -^n 5 " t, l 5 • • • 5 ) 

^det(/;(xf ))^. =1 det( 5l , 2 ^^ 
*^n,3 

Z n , 3 = detA, (43) 
Aij = / fi(x)g 1)2 (x,y)g2,3(y,z)(f)j(z)dfj,(x)dn(y)dn(z) = * g lj2 * 32,3 *4>j- (44) 

JXxXxX 
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Let us write K. and C x as 3 x 3 matrix kernels 

K[x,y)= (45) 

(51,2 * 52,3 * <Pi) ® fj (91,2 * 92,3 * <t>i) ® (fj * 91,2) (91,2 * 92,3 * 4>i) ® (fj * 91,2 * 92,3) 
{92,3 * 4>i) ® fj (92,3 * <Pi) ® (fj * 91,2) {92,3 * 4>i) ® {fj * 91,2 * 92,3) 

4>i ® fj <Pi ® (fj * 91,2) <Pi ® (fj * 91,2 * 92,3) 

51,2 51,2 * 52,3 
- 10 <72, 3 






&(x,y)=J2( A % 1 x (46) 

(51,2 *c 52,3 *c (j>i) <8> fj (51,2 *c 52,3 *c <&) <8> (/j * c 51,2) (51,2 *c 52,3 *c <&) <8> (fj *c 5l,2 *c 52,3) 
(52,3 *c <&) ® /j (52,3 *c <&) ® (/j *c 51,2) (52,3 *c <&) ® (/j *c 51,2 *c 52,3) 

® /j <fo ® (/j *c 51,2) <£i ® (fj *c 51,2 *c 52,3) 

51,2 51,2 *c 52,3 
\ 52 ,3 


Our goal is to show that C x is equal to C x = K-x(Id — K.j)~ l on L 2 (I). The main idea is to use the 
fact that (for arbitrary M ) ICj is "almost" a finite rank operator. Namely, it is equal to a sum of a 
finite rank operator and a nilpotent operator. We first check the identity C 1 = C 1 on a sufficiently 
large finite-dimensional subspace. Let us introduce a finite-dimensional subspace of L 2 (X), 

{( 912 * 523 * (f>i \ I 912 *X 523 * 1 
923 * (Pi , 523 *I 4>i 

V * J V , x „ , . /yj=j <; 

We claim that both /Cj and C 1 leave invariant. 

Lemma 1 The operators K.j, C 1 leave W invariant and C x = C 1 holds on W. 
It follows from straightforward calculations that 




512 * 523 * (p. 

K.J | 523 * (ps 

(ps 

( 912 *1 523 * 4> 
523 *1 (Ps 
o 



( 512 * 523 * (Pi 
523 * (Pi ) X {fj *x 51,2 * 52,3 * 4> S + 

(pi 

' 912 * 923 *1 4>s 






(47) 
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912 *1 523 * 1 
523 *1 4>s 





fj * 91,2 *1 92,3 *1 0s} 



912 * 923 * 9i 

Yl A ij~ I ^23 * H>> 
i,j=l,...n 




x {fj *1 91,2 *1 92,3 *<Ps + 



(48) 




512 *1 923 *1 <Ps 





E V 

i,j=l,...n 

- E A «' 

i,j=X,...n 



912 * 923 * 4>i 
923 * <t>i 



912 * 523 * 4>i 
923 * 4>i 



{fj *X 91,2 * 92,3 *X <Ps}, (49) 
{fj *1 91,2*1 92,3*1 <t>s}- (50) 



Let us introduce the following notations: 



B js = fj *1 91,2 * 92,3 *4>s + fj * 91,2 *1 92,3 * 4>s + fj * 9l,2 * 92,3 *1 
C js = fj *1 91,2 *1 92,3 *(f>s + fj * 91.2 *1 92,3 *1 4>s, 

D 



3S 



E 



fj *1 91,2 * 92,3 *1 4 1 
fj *1 91,2 *1 92,3 *1 • 



j,s= 1, 



,n. 



(51) 
(52) 
(53) 
(54) 



Let us denote by K the matrix of the restriction of the operator /Cj on the subspace W in the basis 

512 * 523 * <Pi \ I 912 *1 923 * 4>i\ ( 912 * 923 *1 4>i \ / 512 *1 923 *1 4> 

923 * <Pi A 923 *1 <Pi , , I . / - 1 II 

<Pi J \ o J \ J \ 
Then it follows from (|47«50|) that 



K 



and 



Id-K 



\ 



/ A- X B A 


- X C A- X D A- 1 


-Id 





-Id 





V o - 


-Id 0. 


Id- A- X B 


-A- X C -A- X D 


Id 


Id 


Id 


Id 





Id 






Id 



J 



The matrix (Id — K) can be easily inverted. The following simple lemma holds. 
Lemma 2 Suppose that a matrix K has a 4 x 4 block form 



( K1 


K 2 


K 3 


K 4 \ 


Id 


Id 








Id 





Id 





V o 


Id 





Id. J 
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and the matrix K\ + — K 2 — K 3 is invertible. Then K 1 is equal to 

/ Q Q(K 4 -K 2 ) -QK 3 -QK 4 \ 

-Q Qi^-Ks) QK 3 QK 4 

-Q Q(K 2 -K 4 ) Ial + QK 3 QK 4 

V Q Q(K 3 - Ki) -QK 3 Id - QK A , ) 

where Q = (K ± + K 4 - K 2 - K 3 )-\ 

In particular for K(Id — K)^ 1 = (Id — K)^ 1 — Id one has the following formula: 



K(Id-K)' 1 



( (A C )~ 1 (A-A C ) -(A c )~ l (E-C) (A c )- l D (A c )' l E \ 

-(A c )~ l A {A c )- l {E-C) -{A c )- l D -(A^E 

-(A^)~ l A {A^)- l {E-C) -{A^)- l D -(A^E 

\ {A )- 1 A — (A°)~ 1 (E — C) — Id {A c )- l D {A c )- l E J 

In the calculations above we used the identity (K1+K4 — K 2 — K 3 ) _1 = (Id— A~ 1 (B+E—C — D)) 
(Id- A^(A- A c ))~ l = (A c )~ l A. 

(91,2 *c92,3 *c < 
92 3 *c 4>i 

can be written as 

512 * 523 * </> 
923 * (pi 

(pi 



912 *X 923 * 4>i 
923 *1 (Pi 





Then 



512 * 523 * 4>s 

C 1 \ g 23 * (p s 

(ps 

912 * 923 * 4>i 
923 * 4>i 



i,j=l,...n 

912 *T 923 * <P: 
923 *1 (Pi 





X {fj *T 91,2 * 92,3 *4>s + fj *c 91,2 *T 92,3 * (ps + fj *c 9l,2 *c 92,3 *T 4>s} 




912 *1 923 * <Ps 
923 *X 4>s 


512 *I 523 * 
& I 523 *1 (Ps 


512 * 523 * <A. 
523 * 4>i 



912 *X 923 *X 4>s 

X {fj *X 91,2 *1 92,3 *(ps + fj *c 91,2 *J 92,3 *1 4>s] ~ ( 
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512 * 523 *J 9. 

C 1 I 



(/ 512 * 523 * <l>i \ ( 512 *1 523 *<Pi \ ( 912* 523 *1 <t>i \ ( 512 *1 523 *1 ' 
523 * 4>i " 523 *X <Pi " + 

\ (pi J \ o / \ o / V o 

{fj *1 51,2 * 52,3 *l4>s}, 
512 *I 523 *1 (Ps 

£ J ( 


(/ 512 * 523 * (Pi \ ( 512 *1 523 * 4>i \ ( 512 * 523 *1 (pi \ ( 512 *J 523 *1 ' 
523 * " 523 *X (Pi " + 

\ <Pi J \ / V / V 

x {fj *X5l,2 *X52,3 *I(ps}- 

It follows from (I5T1-I521) and (JSS1-EH) that £ J = £ z on Lemma is proven. 

To finish the proof of Theorem for the case of three classes of particles, we need to prove the 
relation C 1 = C 1 also on a complement subspace of W in L 2 (l). Let us introduce the following 
subspaces V x C L 2 {h), V 2 C L 2 (/ 2 ), V 3 C L 2 (/ 3 ) : 



Vi 

V-2 

v 3 

We already showed that 



Span{g^ 2 *j 52,3 *i 4>i, i = l,...,n}, 
Span{g 2 ,3 *i 4>i, i = l,...,n}, 
Span{(pi, i = 1, . . . , n}. 



51,2 *X 52,3 *X 9i \ / 51,2 *X 52,3 *X ' 

-J ' 



r J ( I =£ z I I . n. (Gili 



Below we will also prove that 



and 



£ X I 52,3 *X (pi = & 52,3 *i 0i I , i = l,...,n, (61) 

/ V 

C 1 [ I = [ ] . / I n. (02) 



Once this is accomplished, it will be enough to prove C 1 = C 1 on the subspaces 

{Vl)± \ ( ° \ ° \ 

' (^) ± , and of the Hilbert space L 2 (2). What is more, the 

/ V / \ / 

inveribility of the matrices A , A 1,1 , A 1,2 implies that it will be enough to prove the desired relation 
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on the subspaces 



(V4) X \ 


( 







o / 


\ o 



and 



, where V4 = Span{fj,j = 1, . . . , n} C 



L 2 (ii), V 5 = Span{fj *x g lt2 , j = l,...,n} C L 2 (/2) and F 6 = Span{fj *j g lj2 *i 52,3, J = 
1, . . . , n} C L 2 (I 3 ). Indeed, ^ is the matrix of the scalar products of the basis vectors in V\{ resp. 
V 2 , V 3 ) and the basis vectors of V±{ resp. V5, Vq) in L 2 (/i) ( resp. L 2 (I 2 ), L 2 (/ 3 ))). Therefore, 
invertibility of A 1 implies that the sum of V\ and iV^) 1 - is the whole L 2 (I{) (similarly for V 2 and 
(Yb)^, ^3 an d (Vq) ). We claim the following lemma is true. 



Lemma 3 C 1 = C 1 holds on 




04) x \ 


/ ° 


. 




/ 


V 



and 




and on 



92,3 *x 1 









1, . . . ,n. 



The first part is trivial. Indeed, jC 1 = C 1 = on 




To prove the second part of the 



h 





lemma we note that the last identity together with ()6U|) imply that C 1 = C 1 on all vectors 

£ x I 023 *ri I , i = 1, ■ ■ ■ , n, we obtain 



Since it is also true that C x 


that C 1 I 52,3 *x 4> 





512 *1 523 * 4>i 
923 *X 4>i 





By the argument presented above, in order to prove 



C x = C x on all vectors \ h we need to check this relation on (V5) . Let fj *x9\,2 *xh = 0. 



Then 




»J=1 



(63) 



512 * 523 * 4>i 
923 * 




512 *1 923 * <Pi 
923 *I (Pi 
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To calculate C \ h we write 







Now, 



and 



C x \ h = Kj \ h + C x Ki h (64) 



) \ n I gn * 523 * <f>i \ ( 91,2 *i h 

K x [ h \ = Y,A7 j 1 {f j *g 1 , 2 * I h} I fl23 *«fc - ], (65) 

/ i,j=i V </h / V 



\ n / 512 * 523 * <Pi \ ( 91,2 *T h 

C J 1Cj [ h \ = A^ifj * 51,2 *z h} C x I 523 * ^ I - C 1 i ' ) (66) 

/ 5i,2 *x h 

Since fj*xgi,2*xh = we can claim L x \ =0. Substituting (jSSJ) into (|66|) , and combining 

V o 
/o\ /o 

(jnnHnni) we arrive at £ x /i = £ J /i 

So far we have established that C x = C 1 holds on all vectors of the form I /12 I • We also 

V / 

512 * 523 * 4>i \ ( 512 * 523 * <t>i \ ( 512 * 523 * 

proved C x ( 523 * 4>i ) = C % \ 523 * 4>% ), i = 1, ... (since ( 523 * 4>i | G 

/ ° I 

W). Therefore we conclude that C x = C 1 holds on the vectors , i = 1, . . . , n. The last step 



in the proof is to show that C x = CP- holds on I I . Let fj *j 51 2 *i 52 3 *x h = 0. Then 
& ( ) = E ( AC )^' 1 {^ * C 51,2 * c 52,3 *x h) x (67) 

\ /i / *J=1 

512 * 523 * 4>i \ I 512 *1 523 *<Pi \ I 512 * 523 *X </>t \ / 512 *J 523 *X 
923 *4>i - 523 *X 0i - + 

/ V / V / V 

51,2 *c 52,3 *zh\ I 

- 52,3 *xh ) . 

/ V 
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Now. C 1 \ J =K X ( J +C X Ki I | and 
/i / \ /i / \ /t 

\ n I 512 * 523 * <fo \ / 51,2 * 52,3 *jh\ I 

^ J ( = X) ^ 1 {/i*5l,2*52,3*Z^}x 523 *0i - - 52,3 *X h 

h J m=i V <k J V / V 



It follows from (|68|) and previous calculations that CF = C 1 holds on the vectors K-x I I . There- 



h 



fore 



\ / \ n / 512 * 523 * <Ai 

£ J ( =£x + A${ fj* 51,2* 52,3 x£ T 523 *0i 

/» / V / iJ'=l \ <Ai 

51,2 * 52.3 *xh\ I 

C 1 \ 1 - & I 52,3 *x /i 

After simple algebraical calculations we obtain the desired identity. The details are left to the reader. 
Lemma is proven. This finishes the proof of the Theorem for M = 3. 



4 General M Case 

The plan of the proof is the same as in the special case discussed above. We introduce a 2 M ~ 1 x n- 
dimensional subspace of L 2 (I), which is invariant under /Cj and C 1 . As before we will denote the 
subspace by W. The basis vectors of W can be divided into 2 M_1 different groups, each group 
consisting of n elements. The vectors in each group will be indexed by i = 1, ... ,n. These 2 M ~ l 
different groups can be put in one-to one correspondence with the subsets of {1, 2, 3, . . . , M — 1}. 
The empty set will correspond to vectors denoted by e°, i = 1, . . . , n, where can be thought as 

a M- column e- 0) = (51,2*52,3*- ■ ■ * 9m-i,m * 92,3 * ■ ■ ■ *9m-i,m * <f>i, 53,4 * • • • * 4>i, ■ ■ ■ , 
In the case M = 3 the 2 3_1 x n = 4 x n-dimensional subspace W was introduced in the previous 

/ 912 * 523 * <Pi \ 

section. In particular, for M = 3 we have = I 523 * <fii I • The remaining 2 M ~ 1 — 1 

groups will be indexed by . . . , l r ), where r = 1, 2, 3, . . . , l± > 1, . . . , l r > 1, l\ + • • • + l r < 
M — 1. Each such r-tuple corresponds to an r — element subset {h,h + fa, ■ ■ ■ ,h + " ' + lr} of 
{1, 2, 3, . . . , M — 1}. The corresponding basis vectors of W will be denoted by e ^ 1 ' , "' lr \ The vector 
e { ll >— A) i 00 k s s i m Q a r to ef 1 ^ defined above. It can be again viewed as an M-column, but now only 
the first M — l\ — fa — ■ ■ ■ — l r components are non-zero. The first component is 51,2 * • • • * g^h+i *x 
9h+i,h+2 * • • • * 9h+h,h+h+i *T 9h+l 2 +l,h+h+2 *■■■&■ The difference with the first component of e (0) 
is that after 5^,^+1, 5z 1+ / 2 ,/ 1+ / 2+ i, etc, (altogether in r places) the convolution symbol * has been 
replaced by the convolution symbol *x- The second component of e( ll ''"' lr ' is 52,3 * • • • * 5^+1,^+2 *x 
9h+2,h+3 * ' " " * 9h+h+i,h+h+2 *J 9h+i 2 +2,h+h+2 *■■■ The difference with the second component of 
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e(°) is that after gi 1+ ij 1+2 , g; 1 +i 2 +i i / 1 +/ 2 +2, etc, (altogether in r places) the convolution symbol * has 
been replaced by the convolution symbol *%. Please note that in the second component the places 
where symbols * has been replaced by *j are shifted by 1 in comparison with the first component. 
The third component is g 3A *■■■* g h+2 j 1+3 *i 9h+3,h+i *■■■* gi 2+2 ,i 2+3 *i - • fa- Again, the places 
where the convolution symbols *x appear have been shifted by 1 in comparison with the second 
component. In a similar fashion we construct the first M — l\ — l 2 — ■ ■ ■ — l r components. All these 
components have exactlly r convolution symbols *j (in the (M — 1\ — l 2 — ■ ■ ■ — / r )-th component the 
last convolution symbol *j appeares in front of fa). We put the remaining l\ + l 2 ■ ■ ■ + l r components 
to be zero (please note that they can not be constructed according to the principle described above 
if we want to keep the number of convolution symbols *x equal to r). 

Since the constructions of such nature could be better understood after playing with some exam- 
ples, we discuss two of them below. 

Example 

e i = (51,2 *X 92,3 * 93,4 ■■■ * fa, 92,3 *I 93,4 * 94,5 ' " * fa, 9U-1,M *X fa, °)*- Tfle last 

component of is zero (and not, say, fa) because we insist that all non-zero components must 
have the same number of *x convolutions (in this example the number of *x convolutions is one; the 
place of the *z convolution is shifted to the right by a unit each time we go from k-th to k + 1-th 
component, k = 1, 2, . . . , M — 2.) 
Example 

In the case M = 3 (see section 2) the basis vectors ofW are 



912 * 923 *fa\ ( 912 *l 923 * fa 

*l 




ef= 923 * fa , e«= 923 * T fa \ ^ 




912 *X 923 *1 






=1 U I , 1< * < n. 



We start by calculating the matrix of the restriction of ICj on W. It follows immediately from 
23 that 



i=l,...,n l>0 



where 



B (h,...,l r ) _ ^ gi 2 ^ . . . ^ g h h+1 *j g h+1M+2 * ■ • • fa + (69) 

fi * 91,2 *X 92,3 * ' • • 9h+l,h+2 *X 9h+2,h+3 * ' ' ' fa + 

fi * 91,2 * 92,3 *I 93,4 * ■■■ 9h+2,h+3 *X 9h+3,h+A * ' " " fa + ■ ■ ■ , 

where in the first term of the r.h.s. of (|B*9")) the convolution symbols *z appear after fi, gi lt i 2 , gi 1+ i 2j i 1+ i 2+ i, 
(the other convolution symbols are *), in the second term the convolution symbols *j appear after 

9i,2, 9h+i,h+2, 9h+h+i,h+h+2, ■■■ , etc. Altogether, there are M - h - l 2 l r (and not 

M ) terms in the sum, because we require (as before) that each term has the same number of *j 
convolutions. For I + 1\ + • • • + l r > M we agree to set e { s' h '- M = 0. The same rules apply to similar 
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notations introduced below. Let the kernel C x be denned by the right hand side of (|28jl. Then 

£x e {h,-,ir) = ((>l c )- 1 B( ll '-W'% e< (0) ' c -^et I ' ,l --^' c J (70) 

i=l,...,n l>0 



where e ^ 1 ''"'' r ^ c is defined after the example below and q^'—M' i s defined in a similar way to 
g(h,...,l r ) ^ w j^]j a twist. Namely, 

B^ 1 '"'' 1 ^' = fi *i 51,2 * • • • 9h, h+i *i 9h+i,h+2 *---*4>j (71) 

+ fi *c 91,2 *1 92,3 * ■ ■ ■ ^1+1,^1+2 *Z ffli+2,«i+3 * - ' ' 4>j 

+ /j *c 91,2 *c 92,3 *I 53,4 * • • • 9h+2,h+3 *X 9h+3,h+A * ' ' ' 4>j 

+ .... 

The first term of the sum Q71JI is the same as the first term of the sum ()69|) , The only difference 
between the second term in (|69[1 and the second term in (|71|) is that in the second term of (|71|) the 
first convolution symbol (between fi and 51,2) is * c , and not *. In the third term of (|71[) the first 
two convolution symbols are * c and the other are the same as in the third term of (|69|) . etc. 
Example 

Verify that B^< c = A- A c . 

To make sense of (|70|l we also have to define e- ^, and, in general, e^ lv "' ir ^' c . We write 

e i 0) ' C = (51,2 *c • • • *c 9M-\,M *c 4>i,92,3 *c ■ ■ ■ *c 9M-1,M *c <f>i, ■ ■ ■ , <f>iY ■ (72) 

In other words the difference between e[°^' c and e|°^ is that in ef^' c all convolution symbols * are 
replaced by the convolution symbols * c . To obtain e( l i'—M> c from e(h>—M we have to replace in 
each component of e ('i.- -,'r),c first (from the left) l\ — 1 convolution symbols * by * c (in other 
words we do it untill we meet the first symbol *j, at which point we stop). 
The inclusion-exclusion principle implies: 



e (0), c = e (0) + £ £ (-1)^.-^, (73) 



r=l,2,... Zi,...,Z r 

where the summation is over all possible 1 < Zi, . . . , l r , l\ + • • • + l r < M — 1. Similarly, 

e (Zi,...,Z r ),c = e (h,...,J P ) + ^ ^ (_iy e (ki,...,k p ,t,h,-,lr) ^ (j^ 

p=l,2,...,/i — 1 A;i,...,fcp 

where t = Zi — k\ — • • • — k p , and the summation is defined over all possible k\, . . . , k p , such that 
l<h,...,k p , h + --- + k p <h.ln particular, ef u ~ Jr) ' c = e f l '"' M for h = 1. 

Lemma 4 The operators JCj, C x leave W invariant and = C 1 holds on W. 

We have to show that {Id + C T )(Id — ICj) = Id on W. By linearity it is enough to check the 
identity on the basis vectors. It follows from (|68[) and (|7U|) that 



(Id - /Cx)^ 1 -"^ = - (A^B^'-'^iseP + e& h '-' 1 ^ + e (<i>-W (75) 

i=l,. ...n />0 
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and 



(Id + &){Id - K x )e^'-^ = - Y,( A ~ l B {h, - M )iszT ] + E 4 Ml '-'' r) + ei' 1 ""'^ 

i=l l>0 



0).c 



3=1 JS t=l p>l 



+ EE ((A )- 1 ^^'-'^) ef' c - J] e 

3=1 Z>0 ■ ?S p>0,Z>0 

n 



(p,Z,Zl,...,i r ),c 



+ ((A c )- 1 Bf 1 --W' c ) . ef' c _^ e Wi.--M,c (76) 



3=1 JS p>0 

Using the identity B^> c = A - A c we can rewrite (JZHJ) as {Id + C J )(Id - Kx)e^'"' M = S 1 + S 2 + S 3 , 
where 

n M-l 

S X = ^(A-'B^U-e? + e ; 0) ' c + £ ep), (77) 

4 = 1 p=l 



It 



S 2 = Y, i-(A c )- 1 B^ 1 '---' l ^ + {A c )- 1 B^---^' c + Y^ (A c )- l B^~^> c \ e} 0) ' c (78) 

3=1 V Z>0 / j S 

5 , 3 = ^e' l, ' 1, '"' lr ' + e ('l,-W _ e (P' l ' ll '-> lr )' c — ^ e Wl>-,ir),c_ 

«>0 p>0,/>0 p>0 

We claim that Si = S 2 = and S3 = e? 1 '"*'^. Indeed, an easy inductive argument (similar to the 
one showing that B^' c = A — A c ) gives 

M-l 

E = ^ - e S 0) ' c > ( 8 °) 
p=i 

^ £j(Ml>-,Jr),C _ g{h,...,l r ) _ g(h,-,lr),C_ (g-Q 
Z>0 



To tackle S3 we note that 



e (l,h,...,l r ) _ e (l,h,...,l r ),c _ ^2 e (p> l -P>h,-,lr),c^ ^g2) 

1<P</-1 

E E e, w_p ' ,1 '"^' c = Yl e . W, ' 1, '" ,,r),c > ( 83 ) 

Z>0 l<p<Z-l p>0,£>0 

which immediately implies S3 = eg 1 '*" . Lemma is proven. 

To finish the proof of the theorem we need to show that C x = C x on a complement subspace 
of W . We proceed the same way as in section 3. Let us introduce the following subspaces V% C 
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L 2 (h), V 2 C L 2 (/ 2 ), . . . , V M C L 2 (I m ) : 

V\ := Span{g 1)2 *x 52,3 *x 53,4 *z • • • *x <t>i, i=l,...,n}, 
V 2 := Span{g 2 ,3 *i 53,4 *x ■ ■ ■ *x fa, i = l,...,n}, 

V M ■= Span{fa, i = l,...,n}. 

Consider a vector e^' 1 '" (i.e. l\ = l 2 = . . . = Im—1 = 1 ) which has the first com- 
ponent g\ :2 *j (72,3 *J 53,4 fa, and the other components zero. We already proved that 
CF'ef'' 1 '" = C^-ef' 1 '"' , since ef~' ly " € W. In order to prove C x = C x on all vectors of the 
form (h, 0, . . . , o) 1 it is enough to prove the relations for h orthogonal (in L 2 (I\) ) to g\, 2 *j 
52,3 *x 53,4 *x ■ ■ ■ *x fa- The invertibility of matrix A x implies that it is enough to prove it for 
h orthogonal to fj, j = l,...,n, which is trivial since both C 1 and C 1 are identically zero on 
such (h, 0, 0, . . . , 0)*. We proceed now by induction. Suppose that we have already established 
C 1 = C x on the subspace L 2 {h) • • • L 2 (/ m _i) 0{O} • • ■ {0} ( i.e. on the vectors of the form 
(hi, h 2 , . . . , h m -i, 0, . . . ,0)*), where 2 < m < M. We will deduce that the same identity holds 
on L 2 (/i)0---0L 2 (/ m )0{O}---{O}. Consider the vector e ( m - 1 ' 1 >-' 1 ) ) ( i. e . i x = m> ^ = 
... = lu-m-i = 1 )• Since this vector belongs to W, we have /^e^"' 1 ' 1 '-' 1 ) = £l e ( m i 1 > 1 >—' 1 ) < 
The m-th component of this vector is equal to g m ,m+i *x ■■■ *x fa- The inductive assumption 
then implies that C 1 = C 1 on the vector (0, 0, ... , g m ,m+i *x ■ ■ ■ *x fa, • ■ ■ , 0)*, i = 1, . . . , n. 
Using the invertibility of the matrix A 1 and arguing as above, we obtain that in order to prove 
C 1 = C 1 on {0} • • • L 2 (I m ) • • • {0} it is enough to establish the relation only on the vectors 
(0, . . . , h, . . . , 0)*, where only m-th component is non-zero, and fa *j gi )2 *x • ■ ■ *x g m -i,m *x h = 0. 
We have 

/Cx (0,...,/i,...,0)* = Yl A 7 j 1 {fj*9i,2*---*g m -i, m *xh}e[ 0) 

«j=l,...,n 

- (9i,m*xh, g2,m*xh, ...,g m ^ 1;m *xh,0...,Q) t , (84) 
where we remind, the reader that Qi^m — 9i,i-\-i ** * ■ * Sm—i,m 

for 1 < I < m < M. It follows from (|8i|) 
and the inductive assumption that C 1 (/Cx(0, . . . , h, . . . , 0)*) = C x (/Cj(0, . . . , h, . . . , 0)*) . Therefore 

£ J (0,...,/v..,0)* =£r(0,...,fc,...,0)* + & (lCx(0, . . . , h, . . . , 0)*) (85) 
= ^ A i 3 l {fj * 51,2 * • • • * 5m-l,m *J ^}e- 0) - (0i, m *z h, 5 2 , m *Z • • • , 9m-l,m *X h, . . . , 0)* 



+ {{A c )- 1 {A-A c ) A- l ) k . {f j *gi,2*---*gm-x,m*xh} 

k,j=l 



(0),c 



i,j=l,...,n \l>l 



(l),c 



- Y (A^i/ifj *i 9i,m *x h + fj * c 51,2 *x 92,m *x h + fj * c g 1:2 * c g 2j3 *j g^ m *x h + . . .} x e\ 



,(0),c 

(86) 
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+ (51,2 *1 92,m *I h, 0, . . . , of + (g-1,2 *c 52,3 *1 93,m *1 h, 52,3 *1 53,m *lh,Q,..., Of + 

+ (#1,2 *c 52,3 *c 53,4 *J 54,m *J h, 52,3 *c 53,4 *I 54,m *1 h, 53,4 *I 54,m *I h, 0, . . . ,0f + . . . 

+ (51,2 *c 52,3 *c • • • *c 9m-2,m-l *1 9m-l,m *T h, ... , g m -2, m -l *1 9m-l,m *1 h, 0, . . . , 0)* 

(the third term of the r.h.s. of (JHHJ) can be simplified since (A )' 1 (A - A c ) A" 1 = [A c )~ l - A~ v ). 
In a similar way 

& (0, . . . , h, . . . , Of = ( A % l {fj * C 51,2 * c • • • * c g m -i,m *i h}e^' c 

i,j=l,...,n 

~ (91,2 *c 52,3 *c ' ' ' *c 9m—l,m *X h, 52,3 *c ' ' ' *c 5m-l,m *ltl,..., 9m-l,m *X h, . . . , 0)*. (87) 

To see that right hand sides of (|85|) and (|87j) coincide we note that the following three identities hold. 
The first one : 

(51,2 *1 92,m *lh, 0, . . . , 0)* + (51,2 *c 52,3 *I 53,m *1 h, 52,3 *1 53,m *1 K 0, . . . , 0)* 
+ (51,2 *c 52,3 *c 53,4 *J 54,m *1 K 52,3 *c 53,4 *1 9A,m *1 k, 53,4 *J 5 4 , m *J h, 0, . . . ,0f + . . . 
+ (51,2 *c 52,3 *c ■ ■ ■ *c ffm-2,m-l *I 9m~\,m *T h, . . . , 3m-2,m-l *J 5m-l,m *J /l, 0, . . . , 0)* 

= (#1,771 *X /i, 52,m *jh,..., g m -i,m *i h, . . . , 0)* 

- (51,2 *c 52,3 *c • • • *c fa-l,m *J h, 52,3 *c ' " " *c 5m-l,m *i h, . . . , g m -i, m *x h, . . . , 0)*. (88) 

The second one : 

fj *1 5l,m *ih + fj * c 51,2 *J 52,m * J h + /j * c 51,2 *c 52,3 *X 53,m *J h + ... (89) 
+ /i *c 51,2 *c • • • *c ffm-2,m~l *I fe-l,m * I h = fj * 5i, m *jh- fj * c 51, 2 * c ' * * *c 5m-l,m *X ^- 

And the last one is Y^I=i e ^' C = e l°^ ~~ e !°^ C ' which we already used before. The second identity 
allows us to cancel the terms which have the coefficients (A c )^ , the third identity allows us to 
cancel the terms which have the coefficients Aj- , the first identity allows us to cancel the terms 
which contain no coefficients of the form ( J 4 C )~ 1 , AZ ■ All three identities have simple inductive 
proofs relying on a telescoping property of the sums in question. We show the proof of (|89[). The 
proofs of the other two identities are very similar. The main (simple) observation is that 

fj*c9l,2*c- ■ ■*c9m-2,m-l*T9m-X,m*lh+ fj*c9l,2*c ' ■*c<7m-l,m*xft = fj*c9l,2*c ' •*c5m-2,m-l*5m-l,m*J^- 

Taking the sum of the r.h.s. of the last formula and the second to last term of the l.h.s. of (|89j) we 
obtain 

fj *c 51,2 *c • • • *c 5m-3,m-2 * 9m~2,m-l * 9m-X,m *T h. 

At the next step we sum the obtained expression with the third to the last term of the l.h.s. of (|89|) . 
etc. Repeating the procedure the appropriate number of times we obtain fj * 51,2 * • • • * g m -i,m *i h, 
which finishes the proof of (J89|) . Theorem is proven. 

We have learned recently from Harnad that he was able to generalize our theorem to dualization 
with respect to measures modified by arbitrary sets of weight functions (|17j). 

Acknowledgements. It is a pleasure to thank the referees for useful suggestions and John 
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